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Observations of the cosmic microwave background (CMB) and large-scale structure (LSS) provide 
a unique opportunity to explore the fundamental properties of the constituents that compose the 
cosmic dark radiation background (CDRB), of which the three standard neutrinos are thought to be 
the dominant component. We report on the first constraint to the CDRB rest-frame sound speed, 
c^ ff , using the most recent CMB and LSS data. Additionally, we report improved constraints to 
the CDRB viscosity parameter, cj is . For a non-interacting species, these parameters both equal 
1/3. Using current data we find that a standard CDRB, composed entirely of three non-interacting 
neutrino species, is ruled out at the 99% confidence level (C.L.) with c^ ff = 0.30lQ Q2g and c^ is — 
0.44^0 21 (95% C.L.). We also discuss how constraints to these parameters from current and future 
observations (such as the Planck satellite) allow us to explore the fundamental properties of any 
anomalous radiative energy density beyond the standard three neutrinos. 



I. INTRODUCTION 

A complete understanding of the basic building blocks 
of the universe hinges on understanding the elusive prop- 
erties of neutrinos. As a result of having extremely weak 
interactions, neutrinos are the least accurately measured 
of the known fundamental particles. However, even with 
our limited knowledge of neutrino properties, the fact 
that they are massive pQ represents one of the most 
significant challenges to the Standard Model of parti- 
cle physics. Therefore, any improved knowledge of the 
properties of neutrinos will not only serve to shed new 
light on a relatively poorly explored aspect of fundamen- 
tal physics but may also provide further evidence of inad- 
equacies of the Standard Model. Adding urgency to the 
exploration of the properties of neutrinos, a recent com- 
bination of data on anomalous neutrino mixing [5] along 
with observations of the neutrino flux from nuclear reac- 
tors [3] has indicated anomalous mixing between neutrino 
flavors. 

An important consequence of these modifications to 
the neutrino sector is their imprint on cosmological ob- 
servations. Within the Standard Model, a cosmologi- 
cal background of neutrinos is in thermodynamic equi- 
librium with the other cosmological species for tempera- 
tures T v > 1 MeV, after which the neutrino background 
decouples and only interacts through gravity. Neutrinos 
are thought to comprise a significant fraction of the ra- 
diative energy density during big bang nucleosynthesis 
(BBN) causing a measurable effect on the abundances 
of the primordial light elements [4] as well as during the 
formation of the cosmic microwave background (CMB) 
and large scale structure (LSS) [5]. 

Cosmological observations are able to place precise 
constraints on the effective number of neutrino species, 
iVeff, defined so that the total radiative energy density is 
given by p I3d = p 7 [l + iVeff (7/8)(4/ll) 4 / 3 ], where Pl is 
the energy density in photons. The cosmological radia- 



tive content in addition to the photons is known as the 
cosmic dark radiation background (CDRB). In the stan- 
dard cosmological model the only radiative energy den- 
sity besides photons are the three known neutrino species 
so that N e s — 3.046, with the small correction due to fi- 
nite temperature QED effects and neutrino flavor mixing 

Greatly adding to the intrigue, a combination of the 
most current observations of the CMB and LSS indicate 
that N eS > 3 at the 99.9% C.L., with N cS = 4.0±g;|| 
(95% C.L., see Table I). Cosmological constraints on N e g 
are insensitive to anything but the total energy density 
contained in the CDRB. Most radiative backgrounds, 
including neutrinos, share the property that they are 
'non-interacting', i.e., they only interact with other cos- 
mological species through gravity. Examples of such 
backgrounds include axions [7] and short-wavelength 
gravitational- waves [5]. Any further interpretation of 
finding N c g > 3 requires information on other proper- 
ties of the anomalous radiative energy density. 

Here we constrain the values of two additional pa- 
rameters which determine the properties of the CDRB: 
the rest- frame sound speed, c e g, and a viscosity param- 
eter, c v is [S]- As we will describe further, a standard, 
non- interacting, radiative background has (c^ ff ,c^ is ) = 
(1/3,1/3). However, if the CDRB is composed of 
any non-standard species with significant interactions 
these parameters can take on different values (see, e.g., 
Refs. |TQ1 HT] ) which can have a significant impact on the 
observed CMB power spectrum, as shown in Fig. [5] If 
both c^ is and c^ s are found to be consistent with their 
standard value of 1/3 this would lend weight to the in- 
terpretation that observations indicate the existence of 
extra relativistic, non-interacting (i.e., neutrino-like) de- 
grees of freedom. On the other hand, if either is found to 
be inconsistent with their standard values, any inferred 
anomalous radiative energy density cannot be composed 
of standard neutrinos (and may actually be the result of 
unaccounted for systematic effects). 
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Several previous studies have looked at observational 
constraints on neutrino interactions |10[ I12[ 113] . In par- 
ticular, Refs. [13] have used measurements of the CMB 
and LSS to constrain the value of c V i S - However, this 
is the first study to place a constraint on the rest-frame 
sound speed, c c g. Constraints on c c s are particularly 
interesting since models of neutrino interactions indi- 
cate it can differ from its canonical value by ~ 30% 
(e.g., Ref. [10]) and current observations can constrain its 
value to ~ 10% at the 95% confidence level (C.L.). Addi- 
tionally, we report significantly improved constraints on 
c v i s by using the most recent measurements of the CMB 
and LSS lowering the uncertainty on c^ is by a factor of 
-1.5. 



II. PARAMETERIZATION 
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The modified evolution equations for the neutrino per- 
turbations 1 are [9] 
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where the dot indicates a derivative with respect to con- 
formal time, a is the scale-factor, k is the wavenumber, 
c G ff is the rest-frame sound-speed, c V i S is a viscosity pa- 
rameter, b v is the neutrino density contrast, q v is the neu- 
trino velocity perturbation, 7r„ is the neutrino anisotropic 
stress, F Vt i are higher order moments of the neutrino dis- 
tribution function, the shear, a = l/(2k)(h + 6fj), h and r\ 
are the scalar metric perturbations in synchronous gauge 
|14j . and the higher order moments of the distribution 
function are truncated with appropriate boundary con- 
ditions (see, e.g., Ref. [15]). 

We note that the modified evolution equations imply 
a modified set of initial conditions for the perturbation 
equations since neutrinos are a significant fraction of the 
total radiative energy density at early times when the ini- 
tial conditions are set. Following the derivation outlined 
in Ref. [2] we set the initial conditions to the growing 
mode which reverts to the standard adiabatic mode when 



u cff 



1/3 as shown in Appendix A. 



Varying c v i s modifies the ability for neutrinos to free- 
stream out of a gravitational potential well [5] H] . When 



1 Although in this section we refer specifically to neutrinos, these 
equations apply without modification to any massless cosmolog- 
ical component. 



FIG. 1: The evolution of the neutrino density perturbation 
in its rest-frame for a mode with k = 0.1 /iMpc -1 where h 
is the Hubble parameter in units of 100 km/(s Mpc), as a 
function of the scale-factor, a, or the conformal time, r. The 
black solid curve gives the evolution for the standard case, 
i.e., when c^ is = cj? ff = 1/3. The left-hand panel shows the 
evolution when cjj s = (red, dot-dashed) and c^, is = 1 (blue, 
dashed) with c^ ff = 1/3. With c^ ia = (red, dot-dashed) the 
CDRB becomes a perfect fluid leading to undamped acoustic 
oscillations. The right-hand panel shows the evolution when 
c eff = 0.1 (red, dot-dashed) and c^ ff = 0.8 (blue, dashed) with 
c vis = 1/3. When c^ ff is small (red, dot-dashed) the CDRB is 
partially able to overcome its internal pressure support and 
nearly cluster. The bottom panels show the corresponding 
evolution of the Newtonian potential, 3>jv. 
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FIG. 2: Modifications to the CMB temperature power- 
spectrum, C ; TT , as both <4 ia (top panel) and Cgg (bottom 
panel) are varied in the same way as in Fig. [I] the black solid 
curve gives the evolution for the standard case; the top panel 
shows Cf T when c^, ia = (red, dot-dashed) and c^ is = 1 
(blue, dashed); the bottom panel shows — 0.2 (red, dot- 
dashed) and Ceft = 0.7 (blue, dashed). The large angular scale 
measurements are from the 7-year WMAP release [T7] and on 
small angular scales from ACT |19| . 



3 



Cy is = the CDRB becomes a perfect fluid and is capa- 
ble of supporting undamped acoustic oscillations, shown 
in red (dot-dashed) on the left-hand panel of Fig. [I] An 
increased c v i s leads to an overdamping of the perturba- 
tions, shown in blue (dashed) in the left-hand panel of 

Changing c e g allows for a neutrino pressure pertur- 
bation which is non-adiabatic, i.e., (Sp — Sp/3)/p — 

(c 2 s — l/3)6u CSt \ where ^ rest ^ is the density perturba- 
tion in a frame where the neutrino velocity perturbation, 
q v = 0. A value of c 2 s < 1/3 (c 2 s > 1/3) leads to a 
decreased (increased) pressure for the CDRB in its rest- 
frame, which in turn causes the amplitude of the neutrino 
density perturbations to increase (decrease), as seen on 
the right hand side of Fig. [T]in the red, dot-dashed (blue, 
dashed) curve. 

This parameterization is related to a scenario in which 
neutrinos have a significant interaction cross-section. As 
an example, if neutrinos tightly couple to a perfect fluid 
then c 2 ia = and an analogy can be made between our 
parameterization and the tightly coupled photon-baryon 
fluid with a constant sound speed, c 2 , related to c 2 s 
through, 3c 2 s a (c 2 cS + 2/3) . 

We show how the CMB temperature power-spectrum 
is modified in this parameterization in Fig. [2j Note that 
an increase (decrease) in c e ff leads to an increase (de- 
crease) in the scales at which the neutrino perturbations 
affect the CMB. This is due to the increase (decrease) in 
the the neutrino sound horizon. These parameters have 
a similar effect on the polarization power-spectrum, not 
shown here. However, since the effects of the CDRB per- 
turbations are negligible by the time large-scale structure 
forms, the change to the matter power-spectrum is neg- 
ligible [9]. 

III. RESULTS 

In order to measure these parameters we used a mod- 
ified version of the publicly available Boltzmann code, 
CAMB [T5] along with the publicly available Monte Carlo 
Markov chain code, CosmoMC [16] . We used a combi- 
nation of CMB and LSS data including WMAP7 [IT], 
ACBAR [HJ, ACT [19], SPT [20], the Sloan Digital Sky 
Survey (SDSS) DR7 LRG matter power spectrum [2"T] . 
the SDSS small-scale matter power-spectrum measured 
from the Lyman-alpha forest [22] an d the latest determi- 
nation of the Hubble parameter, H , using the Hubble 
Space Telescope [25] . 

In addition to various combinations of (N c g , c 2 is , c 2 s ) 
we allowed the standard six cosmological parameters, 
(A a ,n s ,T,9,fli,h 2 ,fldmh 2 ), to vary, where A s is the am- 
plitude of the primordial power-spectrum, n s is the spec- 
tral index, r is the optical depth, 9 is the angular acoustic 
scale of the CMB, h is the Hubble parameter in units of 
100 km/(s Mpc), fif, is the baryon density in units of 
the critical density, and Sl^m is the dark matter density 
in units of the critical density. All constraints, except 



where noted, force the Helium fraction, Y p , to be fixed 
by its BBN relationship to Qbh 2 and N e s [24]. We note 
that this parameterization only takes into account how 
a change in N e g causes a change in the expansion rate 
during BBN. If the change in N e g is due to a change 
in the physics of the neutrino sector the functional form 
Y p (N e ff) may not hold. We also consider a case where Y p 
is an additional free parameter as discussed below. 
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FIG. 3: Two dimensional contours (68% and 95% C.L.) show- 
ing the degeneracy between c v i s /n s and Ceff/JV e ff. The dot- 
ted contours show the constraints when only CMB data is 
used. The blue-dashed contours show the constraints when 
restricting the CMB to just WMAP7 and large-scale struc- 
ture data, excluding the Lyman-alpha data. The thin-solid 
contours show constraints when restricting the CMB to just 
WMAP7 and with all large-scale structure data. The thick 
contours show the constraints when all of the data are in- 
cluded. The blue circle shows that the standard value of 
c cff — 1/3 and N e s = 3 is excluded at slightly higher than 
the 95% C.L. 

A summary of our results can be found in Table [I] 
Varying the number of effective neutrino species, N c g, 
with c 2 ia — c 2 s = 1/3 we recover a preference for an 
anomalous radiative energy density at the ~ 99.9% C.L. 
When we allow c v - ls and c G ff to vary as well, the signifi- 
cance of any anomalous radiative energy density changes 
to 97.2% C.L. Additionally, the marginalized mean value 
of c 2 s = 0.31 ± 0.015 is lower than the expected value of 
1/3 at the 87.5% C.L. Therefore, as shown by the blue 
circle in the right-hand panel of Fig. [3j the standard value 
of (A G g,Cg ff ) = (3, 1/3) is still disfavored at higher than 
the 95% C.L. 

Since a change in c c g introduces a new length-scale 
(the neutrino sound-horizon) its affect on the CMB is 
only slightly correlated with the other parameters. This 
scale-dependence is clearly shown in Fig. [2] around the 
first peak. Because of its lack of strong correlations, the 
observations place a precise constraint on c 2 s as shown in 
Fig. [5J Constraints to c^ is are not as precise because its 
affect on the CMB is scale-free leading to a degeneracy 
with, for example, the scalar spectral index, n s , shown 
on the left-hand panel of Fig. [3] 

Given the tentative evidence for an anomalously large 
radiative energy density, it is of interest to consider the 
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TABLE I: 

Marginalized ID constraints 
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NOTES.— Errors are 68%, 95% C.L.; upper limits are 95% C.L. 

t For these chains, the number of effective relativistic degrees of 
freedom with c 2 is = c 2 ff = 1/3 is fixed to three while c 2 is and c 2 ff 
applied to one or two extra degrees of freedom are allowed to vary. 



case where the number of neutrinos (with c 2 is 
1/3) is fixed to three and to constrain the values of . 
and 
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applied to an additional one or two effective 
species. We show the results of this analysis in Table 
I. 

Since both c v ; s and c e g modulate the amplitude of 
small-scale power in the CMB we explored degeneracies 
with the running of the spectral index, a s = dn s /dhik. 
Excluding the Lyman-alpha data and fixing c 2 is = c 2 s = 
1/3 we find, a s = -0.020 ± 0.013 ± 0.026 and N cS = 
3.53 ± 0.21 ± 0.72. Allowing both c v i s and c e g to vary we 
find a s = -0.019 ± 0.016 ± 0.03, N eS = 3.49 ± 0.20±g;£g, 
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Although there are many radiative backgrounds which 
are unrelated to neutrinos, if the anomalous CDRB is 
explained by a modification to neutrino physics this may 
lead to a change the functional form of Y p (N c g). There- 
fore it is important to consider the case where Y p is an 



additional free parameter. Fixing c 2 



1/3 we 
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find Y p = 0.294 ± 0.033±g;$$ and N cS = 3.64±g:|i±g:|| 
which is in slight (~1.5<r) disagreement with astrophys- 
ical measurements ^Z5\ and BBN predictions for Y p [2~4] . 
However, when we also allow c v - ls and c or f to vary, the 
helium fraction preferred by the data decreases to Y p = 
0.257 ± 0.051 ±0.1 which is in good agreement with the 
most recent astrophysical measurements, Y p = 0.2565 ± 
0.0010 (stat) ± 0.0050 (syst) [23]. We find the number 
of effective neutrino species is N e g = 3.73+0.28-0.89 w ^h 

^2 



n OO+0.04+0.22 j 2 
u -°°-0.06-0.16 dnu L cff 



0.315 ± 0.018±g:g||. Us- 



ing the astrophysical measurement of Yp as a prior we 
find N cS = 3.73 ± 0.2 ± 0.7, A 
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cf ff = 0.313 ±0.014 
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We note that our modified perturbation equations only 
apply to massless degrees of freedom. However, since a 
non-zero mass predominately affects the late-time (post- 
recombination) evolution of the perturbations, its effects 
are separated in time (and hence we expect should be 
fairly uncorrelated) from the effects of varying c v - ls and 
c e ff , which are most important before and during recom- 
bination. We leave a simultaneous constraint on the neu- 
trino mass, c v i s and c e g to future work. 



IV. CONCLUSIONS 

We have used current CMB and LSS data to explore 
various properties of perturbations in the CDRB. In par- 
ticular, we have parameterized the evolution of the neu- 
trino perturbations with two additional parameters: a 
rest-frame sound speed, c 2 s , and a viscosity parame- 
ter c 2 is , which both equal 1/3 for the standard, non- 
interacting, CDRB. With c% is = c 2 s = 1/3 we find 
that current data favors an anomalously large standard 
CDRB energy density at the 99.9% C.L.. When c vis 
and c e ff are allowed to vary the data still shows that 
the standard value (N e g,c 2 s ) = (3, 1/3) is disfavored at 
slightly greater than the 95% C.L. relative to the 2D- 
marginalized contours. 

Our results can be interpreted as providing tentative 
evidence that the extra relativistic degrees of freedom 
seen in observations of the CMB may have non-negligible 
interactions with c 2 s < 1/3 at the 87.5% C.L. Although, 
as shown in Fig. |3j this result is driven by the small- 
scale observations of the CMB and may be impacted by 
systematic errors, such as issues related to marginalizing 
over sources of secondary anisotropics. A more conclusive 
result must wait for data from future observations, such 
as the Planck satellite |26) . 

If any anomalous radiative energy density is due to a 
modification of the neutrino sector this may result in a 
change to the standard relationship Y p (N e g). Allowing 
the helium fraction to also vary we find that the con- 
straint on Y p is in agreement with the most recent astro- 
physical measurements [23J . In addition the values of c 2 s 
and c 2 is are consistent with their non-interacting value 
of 1/3 and iV c ff is larger than the expected value of 3 
at the 90% C.L. Using the astrophysical measurement of 
Y p as a prior we find N e g > 3 at the 95 % C.L., cj is is 
fully consistent with the expected value of 1/3, and c^ s 
is less than 1/3 at the 85% C.L. Therefore we find that 
although fixing the helium fraction through its BBN re- 
lationship Y p (N c g) may not be appropriate in general, 
when we allow Y p to be a free parameter the constraint 
on N c g is still significantly anomalous. Using the astro- 
physical measurement as a prior on Y p only increases this 
significance and hints at a slightly low value for Cg ff . 

Using a Fisher analysis we find that future observations 
of the CMB with the Planck satellite alone will provide 
a measurement of N e g = 3.0 ± 0.17, = 0.333 ± 0.026, 
and Cg ff = 0.333 ± 0.004. If future observations continue 
to provide evidence for the presence of extra relativistic 
energy density then, when applied to one additional ef- 
fective neutrino degree of freedom, Planck will constrain 
4 is = 0.3 ± 0.1 and c 2 cS = 0.333 ± 0.017. The improved 
sensitivity to these parameters will allow a constraint on 
the fundamental properties of any new radiative degrees 
of freedom. 

Note added in proof: after a preprint of this paper 
appeared on the arXiv we became aware of a study, 
Ref. [37], which presents similar results. 
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where ML denotes the massless neutrinos (parameter- 
ized by c e ff and c v i s ), MV denotes the standard mas- 
sive neutrinos, 5i are the density contrasts, qi is the heat 
flux, Hi is the anisotropic stress, z — \h where h is 
the standard synchronous gauge potential (see Ref. [14]). 
R ML = pML /ptot ig the f rac tion of the total radiation 
energy density in massless (non-standard) neutrinos, and 



[1] Y. Fukuda et al. (for the Super-Kamiokande collab- 
oration), P hys. Rev. Lett. 81, 6 (1998) |arXiv:hep-] 
ex/9805021] . 

[2] A. Aguilar et al. (LSND), Ph ys. Rev. D 64, 112007 
(2001) [arXiv: |hep-ex/0104 049 ; A. A. Aguilar-Arevalo 
et al. (MiniBooNE), Phys. Rev. Lett. 105, 181801 (2010) 
[arXiv: 1007.1150 [hep-ex]]. 

[3] Th. A. Mueller et al., Phys. Rev. C 83, 054615 (2011) 
[arXiv: 1101.2663 [hep-ex]]. 

[4] R. H. Cyburt, B. D. Fields, K. A. Olive and E. Skill- 
man, Astropart. Phys. 23, 313 (2005) [arXiv: |astro-| 



Center for Cosmological Physics. 

Appendix A: Initial conditions 

Following the derivation outlined in Ref. [T3] we set 
the initial conditions to the growing mode which reverts 
to the standard adiabatic mode when c 2 is = c 2 ff = 1/3. 
For reference, these initial conditions are given in syn- 
chronous gauge by 

(Al) 
(A2) 
(A3) 

(A4) 

(A5) 
(A6) 

(A7) 

(A8) 
I 

Rf v = p f v /pt°t is the fraction of the total radiation 
energy density in massive, standard, neutrinos. One can 
check that that for standard neutrinos (c 2 ff = c 2 is = 1/3) 
these initial conditions revert back to the standard adi- 
abatic initial conditions; for c 2 ff ^1/3 the initial con- 
ditions are an admixture of adiabatic and isocurvature 
initial conditions. 



ph/0408033 ; K. A. Olive, G. Steigma n and T. P. W alker, 
Phys. Rept. 333, 389 (2000) [arXiv: astro-ph/9905320| . 
[5] W. Hu, D. Scott, N. Sugiy ama and M. J. White , 
Phys. Rev. D 52, 5498 (1995) |arXiv:astro-ph/9505043l ; 
W. Hu and N. Sugiyama , Astrophys. J. 471, 542 (1996) 
arXiv:astro-ph/9510117 ; R. Bowen, S. H. Hansen, 
A. Melchio rri, J. Silk and R. Tro tta, MNRAS 334, 
760 (2002) [arXiv:astro-ph/0110636"l Z. Hou, R. Keisler 
L. Knox, M. Marius and C. Reichardt, [arXiv:1104.2333| 
[astro-ph.CO]]. 

[6] G. Mangano et al., Nucl. Phys. B 729, 221 (2005) 



5c = S b = \&t IV = \5, = -|fcV, 

sML = f 2[10 + <4(2 + 7Rt IL + Rt IV )] - 3c c 2 ff [5 + 2^(2 + Rf L + R$ IV )} , 

10 + 3^ + 6(1 + 3^)4^ < 7 ' 



kr 

It = "g~ T' 



i / 4(2+iff*n 

% ~ 9 v 1+ Ih + AR^v ) krb ' 



j 7 , 



ML ^2c? vis {cl s [6RM L - 3(2 + RMV)}- 2(2 + R*fL + RMV)} 



v v = X 



10 + 3c c 2 ff + 6(l + 3c c 2 ff )4 is i?^ 



mv _ 2(2 + Rf L 2 2 

^ _ X i5 + 4i?^y fc r ' 



ml J 2<4(2 + Rt IL + Rt IV + 3c c V [5 + 2<4(2 + R^ L + R$ IV )\ 
? " 1 30 + 45^ + 18(1 + 3^)^^ f*^ 



6 



arXi v:hep-ph/ 0506164 . [17] 
[7] M. S. Turner, Phys. Rev. Lett. 59, 2489 (1987). 
[8] T. L. Smith, E. Pierpaoli and M. Kamionkowski, Phys. 

Rev. Lett. 97, 021301 (2006) |arXiv:astrc~ph/0603144]. [18] 
[9] W. Hu, As trophys. J. 506, 485 (1998) |arXiv:astro-| 

ph/9801234] . 

[10] N. F. Bell, E. Pier paoli and K. Sigurdson, P hys. Rev. D [19] 
73, 063523 (2006) [arXiv:astro-ph/0511410] . 

[11] R. F. Sawyer, Phys. Rev. D 74, 043527 (2006) [20] 
[arXiv:astro-ph/0601525] . 

[12] S. Hannestad, JCAP 0502, 011 (2005); M. Cirelli [21] 
and A. Stru mia, JCAP 0612, 013 (2006) [arXiv:astro^| 
ph/0607086]; A. Friedla nd, K. M. Zurek and S. Bashin- [22] 
sky. [arXiv:0704.3271| [astro-p h.CO]1; A. F riedland, 
K. M. Zurek and S. Bashinsky. [arXiv:0704.327l1 [astro- [23] 
ph. CO]]; A. Basboll, O. E. Biaelde, S. Hannestad and 
G. G. Raffelt, Phys. Rev. D 79, 043512 (2009) [arXiv: [24] 
0806.1735 [astro-ph]]. 

[13] R. Trotta and A. Melchiorri, Phys. Rev. Lett. 95, 011305 [25] 
(2005) [arXiv:astro~h /0412066 ; F. De Bernardis, 
L. Pagano, P. Serra, A. Melchiorri and A. Cooray, JCAP 
0806, 013 (2008) [arXiv:astro-ph/0804.1925]; 

[14] C. P. Ma an d E. Bertschinger Astrophys. J. 455, 7 [26] 
(1995) [arXiv:' |astro-ph/9506072] . 

[15] A. Lewis, A. Challinor and A. Lasenby, Astrophys. J. [27] 
538, 473 (2000) |arXiv:astro-ph/9911177] . 

[16] A. Lewis and S. Bridle, P hys. Rev. D 66, 103511 (2002) 
|arXiv:astro-ph/0205436] . 



D. Larson et al. (for the WMAP collaboration), Astro- 
phys. J. Suppl. 192, 16 (2011) [arXiv:1001.4635l [astro- 
ph.CO]]. 

C. L. Reichardt et al. (for the ACBAR collaboration), 
Astrophys. J. 694, 1200 (2009) [arXiv:0801.l"49"i"1 [astro- 
ph.CO]]. 

S. Das et al. , Astrophys. J. 729, 62 (2011) 
[arXiv: 1009.0847] [astro-ph.CO]]. 

R. Keisler et al. (for the SPT collaboration), Astrophys. 
J. 743, 28 (2011) jarXiv:1105.3l"82] [astro-ph.CO]]. 
B .A. Reid et a l. (for the SPSS c ollaboration), MNRAS 
404, 60 (2010) [arXiv:09071659| [astro-ph.CO]] . 
P. McDonald et al. (for the SPSS colla boration) , AJ 635, 
761 (2005) [arXiv: [astro-ph/0407377] . 
A. G. Riess et aTJ Xstrophys. J. 730, 



|arXiv:110 3.2976 [astro-ph.CO]]. 
G. Steigman, Ann. Rev. Nucl. Part 
|arXiv:0712.1100| [astro-ph]] . 
Y. I. Izotov, T. X. Thuan and G. Stasinska 
J. 662, 15 (2007) !;arXiv:astro-ph/0702072l 
tov and T. X. Thuan, Astrophys. J. 710, 



119 (2011) 
Sci. 57, 463 (2007) 



Astrophys. 
; Y. I. Izo- 
L67 (2010) 
arXivilOOl 4440| [astro-ph.CO]]; 

3 . A. R. Ade et al. (for the Planck Collaboration), 



arXiv:11 01.2022 [astro-ph.IM]]. 
M. Archidiacono, E. Calabrese, and A. Melchiorri 
jarXiv: 1109.2767] [astro-ph.CO]]. 



